Lagrangian temperature, velocity and local heat flux measurement in 

Rayleigh-Benard convection 
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We have developed a small, neutrally buoyant, wireless temperature sensor. Using a camera for 
optical tracking, we obtain simultaneous measurements of position and temperature of the sensor as 
it is carried along by the flow in Rayleigh-Benard convection, at Ra ~ 10 10 . We report on statistics 
of temperature, velocity, and heat transport in turbulent thermal convection. The motion of the 
sensor particle exhibits dynamics close to that of Lagrangian tracers in hydrodynamic turbulence. 
We also quantify heat transport in plumes, revealing self-similarity and extreme variations from 
plume to plume. 

PACS numbers: 47.80.-v (Instrumentation for fluid flows); 44.27.+g (Convective heat transfer) 



Understanding fluid motion and transport of heat 
due to thermal convection is crucial for progress in di- 
verse challenging and important problems such as climate 
change, processes in planetary and stellar cores, and ef- 
ficient temperature control in buildings. Numerous lab- 
oratory studies employing Rayleigh-Benard experiments 
have uncovered laws which relate global heat flux and 
flow velocities to fluid properties, flow boundary geom- 
etry, and driving parameters [lj. More recently, several 
studies @, [1| have focused on coherent flow structures, 
called plumes, which are ejected from the thermal bound- 
ary layers carrying heat into the convective flow. Wc 
present measurements from a novel temperature sensor 
which is carried along with the convective flow, i.e. La- 
grangian measurements. In this way we obtain statistics 
of simultaneous velocity, temperature, and heat trans- 
port dynamics throughout the lifetimes of many plumes. 

Lagrangian measurements are particularly well suited 
for studying flows where coherent structures Q and 
mixing [5J are important, e.g. plumes mixing temper- 
ature in convection. Recent Lagrangian experiments 
with submillimeter-sized passive tracerparticles have ad- 
vanced understanding of turbulence [a 0, B Q ■ In the 
context of convection, Lagrangian measurements have 
long been performed using meter-sized atmospheric bal- 
loons [3] or ocean floats [Tl[ ■ Here we use a similar strat- 
egy, aided by advances in miniature sensors and com- 
munication devices to probe temperature and flow 
properties at centimeter scales in Rayleigh-Benard con- 
vection. 

Our experimental setup is a rectangular vessel with 
height H = 40 cm and section 40 cm x 10 cm filled with 
water (for more detail, see [HI)- The walls are made 
of poly-methylmetacrylate (PMMA), the top boundary 
is a copper plate chilled by a controlled water bath, 
and the bottom boundary is an electrically heated cop- 
per plate. The heater power is maintained at 230 W 
and the top plate is held at T up = 19°C, resulting in 
a temperature difference AT = 20.3°C between the top 
and bottom plates. The resulting Raylcigh number is 



Ra = g(3ATH 3 /vn = 3x 10 10 
due to gravity, (3 = 2.95 x 10 



where g is acceleration 
4 K" 1 is the thermal ex- 



pansion coefficient of water and v = 8.17 x 10 _7 m 2 s _1 , 
K = 1.48 x 10 m 2 s _1 its viscosity and thermal diffu- 
sivity (values are given for the mean temperature of the 
flow 29.1°C). The Nusselt number, defined as the total 
heat flux normalized by kAT/H, is Nu = 167.9 ± 0.2. 
Under these conditions, the convective regime is fully 
turbulent [U 
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1 1 51 ] and the mean flow is a steady, 
system-sized, convection roll with a rotation period of 
about 100 s. 

The mobile sensor consists of a D = 21 mm diameter 
capsule containing temperature instrumentation, an RF 
emitter, and a battery. It is described in detail in [l2j and 
we recount its basic features here. The capsule and fluid 
density arc carefully matched within 0.05 percent so that 
it reliably follows the flow. Four thermistors (0.8 mm, 
230 kil, response time 0.06 s in water) are mounted in 
the capsule wall protruding 0.5 mm into the surrounding 
flow. A resistance controlled oscillator is used to create 
a square wave whose frequency depends on the tempera- 
ture of the thermistors. This square wave is used directly 
to modulate the amplitude of the radio wave generated 
by the RF emitter. The temperature signal is recov- 
ered on-the-fly by a stationary receiver and a Labview 
program. The dynamic range of temperature detection 
is 80dB, with a resolution of 4 mK and 50 ms. In ad- 
dition, the capsule trajectory is recorded with a digital 
video camera, providing synchronous measurements of 
the position and temperature of the sensor as it is car- 
ried about by the fluid. With maximum flow velocities in 
the range 1-2 cm/s and a particle size of 21 mm, we are 
oversampling the dynamics by a factor of order 10. How- 
ever, the characteristic thickness of the thermal bound- 
ary layer is £t ~ ^HNu^ 1 ~ 1.2 mm and that of the 
hydrodynamic boundary layer may be estimated [l[ as 
lu ~ It(v I k)+ 1 I z ~ 2 mm. Thus, the sensor is too large 
to penetrate the boundary layers. 

Before presenting more detailed statistics, we provide 
an overview of the raw data collected by the mobile sen- 
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FIG. 1: Temperature and velocity measurements of the mobile sensor, (a) trajectory with temperature color coding; (b) 
temperature; (c) vertical position; (d) vertical velocity; (e) cross correlation of temperature and velocity. 



sor. Predominantly, the sensor moves in looping trajecto- 
ries the size of the convection cell with a period of about 
100 s. The collective result of many fluid parcels (plumes) 
with such trajectories compose the large scale convection 
roll - one may also detect in Fig. la the presence of sec- 
ondary rolls in the lower left and upper right corners. 
Although the main motion of the sensor is rather peri- 
odic, we point out that its temperature fluctuates widely 
and irregularly (see Fig. lb, c, d.) Fig.le shows the cross- 
correlation between temperature and normalized veloc- 
ity. We find a maximum value of 0.6 at zero-time lag 
and confirmation of the 100 s period of the large scale 
roll. This correlation is about twice the value reported 
in for local Eulerian measurements. We attribute 
these dynamics to the entrainment of the sensor by ther- 
mal plumes and associate the fluctuations in Figsl.b-d 
with variation between plumes. In addition to qualita- 
tive verification of this idea using Schlicren visualizations, 
it is consistent with previous studies which find plumes 
predominately near the side walls like the trajectory of 
our sensor 16|, [l7[ . We explore these dynamics quantita- 



tively in the results to follow. 

From the data in Fig.l, one may estimate character- 
istic dimensionless numbers for the flow. From inte- 
gral quantities one computes the integral Reynolds num- 
ber Re = UH/v ~ 4000. The associated turbulent 
- Taylor based - Reynolds number is of the order of 
R\ oc \f~Re ~ 63. We note that this value is very close to 
a local Reynolds number, defined from the actual motion 
of the Lagrangian sensor Re 1 = u rms t rms jv ~ 60, where 
l rm s is the fluctuation of position along a 'mean' trajec- 
tory, u rms is the usual standard deviation of the velocity. 
In this regime (R\ ~ 60), one expects turbulence with 
significant intermittency in the velocity gradients 

Accordingly, we turn now to more detailed statistics 



of our measurements including power spectra, increment 
probability distribution functions (PDF), and structure 
function scaling for velocity and temperature. The power 
spectra (in time) are shown in Fig. 2a. The velocity spec- 
trum is close to the behavior expected for a Lagrangian 
tracer in a turbulent flow - f~ 2 - although scaling range 
is limited at such low R\ . The temperature spectra is ob- 
served to roughly mimic that of velocity with a steeper 
power law slope compared to measurements with a sta- 
tionary temperature probe. For such Eulerian measure- 
ments, experiments [19| measure slopes between -1.35 
and -1.4 (= —7/5) while numerical studies [2(| report 
-7/5 in the bulk and -5/3 near the side walls. Veloc- 
ity increments (Fig. 2b) indicate that the sensor is sub- 
jected to rather intermittent acceleration with strongly 
non-Gaussian statistics - flatness is -21. Like the power 
spectrum, the shape of PDFs for different velocity in- 
crements is similar to those of a fluid particle in a tur- 
bulent flow - Fig. 2b. Finally, we have computed the 
evolution of the structure functions for the vertical ve- 
locity S% z (t) = (\v z (t + t) — v z (t)\ p ) t ; we observe an 
extended region of relative scaling S$ z (t) oc [5^(r)] ?p 
(for 0.3 s < r < 10 s) with exponents £i = 0.56,^2 = 
1.00, & = 1.30, & = 1.50, & = 1.65,6s = 1.7. These 
values (determined with a 10% precision) are again in 
good agreement with experimental measurements for La- 
grangian tracers in turbulence 0. We also show in 
Fig. 2c the probability density function (PDF) of temper- 
ature increments measured by the mobile sensor. Similar 



to Eulerian measurements 17[, their statistics are non 
Gaussian, with wider tails at small scale. In contrast to 
the velocity behavior (and to some Eulerian studies [iol]). 
our Lagrangian temperature increments do not reveal a 
range of self-similar scaling. 



As proposed in several previous studies [16|, [2l|, |22|, 
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FIG. 2: Temperature and velocity statistics, (a) Time spectra 
of temperature and vertical velocity of the sensor particle; the 
frequency is non-dimensionalized using the period Tr — 100 s 
of the large scale roll motion and curves have been shifted 
vertically for clarity). The solid straight line corresponds to 
an ideal f~ 2 scaling; (b) PDFs of vertical velocity increments 
(time lags are given in the inset); (c) PDFs of temperature 
increments for time lags equal to 0.1 and 102.4 s. 



the dynamics of heat transport may be best analyzed in 
terms of a local heat flux v z 9'(t), where 9'(t) = 9{t) — 9 
is the particle temperature variation from its time aver- 
aged temperature 8, and v z (t) its vertical velocity. Our 
mobile sensor allows for a Lagrangian measurement of 
this quantity. In particular, we define a normalized La- 
grangian vertical heat transport 



Nu L {t) = 1 + 



H 
kAT 



v z (t)9\t) 



(1) 



The time series Nu L (t) is shown in Fig. 3a. As expected, 
it is most often positive since convective motions are as- 
sociated with either hot fluid rising or cold fluid sinking, 
in each case (v z 9') > 0. The much less probable events 
with (v z 9') < correspond, for instance, to the rise of 
the particle when it is colder than its environment due 
to turbulent swirls in the flow. We find that the time 
averaged Lagrangian heat transport Nu L = 328 is larger 
than the global value (Nu = 168 in our case). This 
indicates that our sensor preferentially samples the re- 
gions of the flow which carry higher than average heat, 
i.e. plumes. In contrast, the traditional Nussclt num- 
ber indiscriminately accounts for all regions in the flow 
such as the center of the cell where plumes rarely visit. 
Another prominent feature is the highly non-Gaussian, 
intermittent fluctuations of Lagrangian heat transport, 
which is similar to measurements with stationary probes 
at similar Rayleigh numbers 17 1. 

To further quantify Lagrangian heat transport, we 



FIG. 3: Lagrangian heat transport, (a) time series; (b) PDF 
of increments, and Gauss distribution (dashed line), (c) ESS 
plot of fourth and sixth order structure functions, vs the sec- 
ond order one (the second order structure function has been 
shifted vertically). 



have studied the statistics of the increments (in time) 
Nu L (t + t) — Nu L (t). Two features are noteworthy: (i) 
as shown in Fig. 3b, their PDFs arc strongly non Gaus- 
sian from the smallest time increments to lags of the 
order of integral time Tr; (ii) there exists a range of 
relative scaling for their structure functions S p Nu (t) = 
(\Nu L (t + t) - Nu L (t)\P) t , as shown in Fig.3c. The 
relative exponents, S^ u (t) oc [^jvuC 7 ")] have values 
£f" = 0.54, ^ = LOO, Cf" = 1.37, Cf u = 1.68, £f u = 
1.93, £e = 2.14 - the use of the second order structure 
function as a reference is arbitrary. Comparison with the 
velocity shows that the Lagrangian heat flux is less inter- 
mittent; i.e. the PDFs of the local heat flux in Fig. 3b are 
non Gaussian but their shapes do not evolve significantly 
as increments increase towards the large scale period Tr. 

A simple interpretation of the above statistical results 
is that thermal plumes may be defined as Lagrangian 
heat transport events (once detached from the thermal 
boundary layer); these events are in some way self-similar 
as will become clear below. In this light, we have ana- 
lyzed the portions of the sensor trajectories during which 
its vertical position lies between if/4 and 3i//4, i.e. the 
rise (or fall) of a hot (cold) plume. For each of these mo- 
tions, we compute the mean heat transport, and its fluc- 
tuation measured by the standard deviation computed 
over the portion of the trajectory. The results are re- 
ported in Fig. 4. We observe that the heat transport has 
large variations from one plume to the next. The mean 
Lagrangian heat transport in these events (Nu pl = 335) 
is quite close to the mean for the entire time series 
(Nu L = 328), which indicates that the times the sensor is 
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not in a plume arc relatively unimportant to heat trans- 
port. Furthermore, the variation from plume to plume is 
quite large, rms(ATu pl ) = 384, of the order as the mean. 
During the trajectory of a single plume, the standard de- 
viation a of Nu pl is proportional to the mean; a = aNu pl 
where a = 0.5 ± 0.05. We note that the plume heat 
transport we measure matches within a factor of 2 the 
prediction of the model proposed in [22| . though more 
experiments are required to test scaling laws for different 
Rayleigh numbers and fluid properties. Furthermore, we 
find that the characteristic time for temperature decay 
during the life-time of a plume is of the order of its travel 
time from one plate to the other, which also supports the 
model in El. 
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useful discussions with Sergio Ciliberto and Romain Volk. 
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FIG. 4: Plumes. Mean (dots) and standard deviation (error 
bars) of Lagrangian heat transport of consecutive plumes. 

To summarize, we have reported novel measurements 
of Lagrangian temperature and heat transport using a 
wireless, neutrally buoyant temperature sensor. The 
sensor provides a new perspective on the dynamics of 
thermal plumes in turbulent Rayleigh-Bcnard convec- 
tion. We find that heat transport fluctuates greatly 
from one plume to another and that these fluctuations 
suggest a self-similar character of plumes. A practical 
implication of our results is that in order to maximize 
heat transport, one should maximize the production of 
thermal plumes. This is supported by global Nu versus 
Ra measurements in experiments 23j, |24J with rough 
cndplates. Future investigations will be focused on the 
scaling of plume heat transport for varying Ra and 
different fluids. 
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